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FIGURES  1-36 
APPENDIX  A: 

CORRELATION  BETWEEN  WAVE  HEIGHTS  AND  ES¬ 
TIMATES  BY  SEISMOMETER  . 
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Part  1 

Introduction 

The  extreme  wave  climate  of  the  Pacific  Northwest  and  the  need  for  nearshore 
wave  data  led  to  an  investigation  into  the  performance  of  the  seismome¬ 
ter  wave  gage  by  the  US  Army  Engineer  Waterways  Experiment  Station’s 
Coastal  Engineering  Research  Center  (CERC),  Vicksburg,  Mississippi,  to 
evaluate  its  potential  for  satisfying  wave  data  requirements  of  the  US  Army 
Engineer  District  (USAED),  Portland,  Oregon.  Data  were  recorded  simul¬ 
taneously  from  the  seismometer  and  from  the  output  of  the  attached  Oregon 
State  University  “wave  meter”  device  during  the  summer  of  1985  at  Brook¬ 
ings,  Oregon.  Additional  data  were  obtained  for  the  same  time  periods 
from  a  directional  wave  gage  installed  by  Oregon  State  near  Brookings  as 
well  as  other  available  data.  The  data  were  analyzed  and  compared,  and 
various  means  of  improving  the  analysis  methodology  were  investigated. 

Wave-height  estimates  from  seismometer  data  can  be  shown  to  have  a 
high  correlation  with  actual  sea  states.  Overall,  however,  the  estimates 
by  the  seismometer  are  less  accurate  than  required  for  wave  statistics  for 
US  Army  Corps  of  Engineers  (Corps)  project  designs.  Attempts  to  dif¬ 
ferentiate  between  local,  nearshore  wave  conditions  and  offshore  or  distant 
conditions  were  partially  successful,  but  additional  comparison  data  and 
work  are  needed.  Improvements  in  wave  period  estimates  were  success¬ 
fully  achieved,  and  reliable  estimates  of  offshore  dominant  wave  periods  are 
probably  obtainable. 

A  seismometer- based  wave  gage  could  be  useful  for  wave  estimates  for 
operational  purposes  where  the  requirements  for  real-time,  low-cost,  reli¬ 
able  data  would  offset  the  limitations  on  system  accuracy.  Wave  data  from 
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such  a  gage  would  infrequently  be  unavailable  due  to  adverse  seismic  ac¬ 
tivity,  and  occasionally  have  large  errors  in  estimates.  There  is,  however, 
some  indication  that  reliability  during  extreme  wave  events,  which  are  of 
most  interest,  would  be  better  than  average.  A  seismometer  coupled  to  a 
local  microcomputer-based  system  could  provide  estimates  via  telephone, 
radio,  or  satellite  links  to  Corps  operational  personnel  or  contractors.  De¬ 
velopment  of  such  a  system  is  well  within  the  state  of  the  art  and  would 
make  use  of  the  analysis  techniques  developed  by  this  study. 


1.1  Objectives 

The  USAED,  Portland,  Oregon,  requires  synoptic  measurements  of  waves 
in  the  coastal  and  nearshore  waters  along  the  coast  of  Oregon  to  perform 
its  various  missions  in  planning,  engineering,  and  operations.  In  the  mid 
1970’s  a  system  for  sea-state  measurement  based  on  the  principle  of  micro¬ 
seisms  was  developed  by  Oregon  State  University  (OSU)  and  deployed  at 
numerous  Oregon  coastal  sites  by  the  National  Oceanic  and  Atmospheric 
Administration. 

The  CERC  has  undertaken  an  evaluation  of  the  existing  OSU  seismome¬ 
ter  wave  meter  to  evaluate  its  accuracy  and  to  determine  what  improve¬ 
ments,  if  any,  could  be  made.  In  the  first  phase  of  this  work,  Thompson, 
Howell,  and  Smith  (19S5)  reviewed  the  design  of  the  wave  meter  system 
and  analyzed  the  data  manually  using  the  wave-by-wave  analysis  method 
of  CERC  to  compare  with  other  wave  data.  The  results  for  several  winter 
storm  periods  indicated  surprisingly  good  agreement  for  significant  wave- 
height  estimates,  while  period  estimates  often  showed  substantial  disagree¬ 
ment.  Several  recommendations  for  improvement  to  the  design  of  the  wave 
meter  were  also  made. 

The  objectives  of  the  present  work  were  to  continue  the  evaluation  of 
the  existing  wave  meter  system  during  summer  wave  conditions  using  digi¬ 
tally  recorded  data  and  computer  analysis  and  to  investigate  microseismic 
analysis  techniques  which  might  be  useful  to  improve  the  seismic  wave- 
measurement  system. 
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1.2  Background 

The  existence  of  small  amplitude  seismic  waves  of  frequencies  between  0.04 
and  0.4  Hz  has  long  been  known.  They  have  been  detected  by  numerous 
observations  throughout  the  world,  mainly  in  coastal  areas.  Attempts  to 
explain  the  mechanism  of  the  high-frequency  seismic  noise  have  appeared 
in  a  number  of  publications.  For  general  review,  reference  should  be  made 
to  Darbyshire  (1962)  and  Hasselmann  (1963). 

The  theory  that  ocean  surface  gravity  waves  can  cause  seismic  waves 
is  well  established  analytically  and  supported  by  many  experiments.  It  is 
well  known  that  standing  surface  waves  lead  to  the  fluctuations  of  the  mean 
dynamic  pressure  on  the  sea  bottom.  In  this  case,  the  pressure  field  has  a 
component  unattenuated  by  water  depth,  and  the  frequency  is  double  that 
of  the  standing  surface  wave.  This  nonlinear  (with  respect  to  the  wave  am¬ 
plitudes)  theory  was  reviewed  and  applied  to  the  problem  of  microseisms 
by  Longuet-Higgins  and  Ursell  (1948)  and  again  later  by  Longuet-Higgins 
(1950)  with  a  complete  microseism  generation  theory.  Cooper  and  Longuet- 
Higgins  (1951)  reported  a  wave  tank  experiment  on  the  pressure  variations 
due  to  standing  waves.  Although  Longuet-Higgins’  explanation  of  micro¬ 
seisms  is  based  on  the  precondition  of  standing  waves  and  has  not  always 
been  supported  by  experiments  (see  Donn(1952)),  many  observations  have 
confirmed  the  relationship  between  the  frequencies  of  ocean  waves  and  mi¬ 
croseisms  and  thus  supported  this  theory  (see  Deacon  (1947),  Dinger  and 
Fisher  (1955),  Carder  (1955),  and  Haubrich,  Munk  and  Snodgrass  (1963)). 

Zopf,  Creech  and  Quinn  (1976)  suggested  a  method  for  ocean  wave  mea¬ 
surements  based  on  seismic  recordings,  employing  the  generation  theory  of 
Longuet-Higgins  (1950).  The  method  assumed  that  relationships  between 
amplitudes  of  the  elastic  deformation  of  the  sea  bottom  and  the  pressure 
caused  by  ocean  waves  are  linear  and  thus  the  amplitudes  of  ocean  waves 
can  be  measured  by  an  empirical  calibration  curve  representing  a  relation¬ 
ship  between  amplitudes  of  microseisms  and  ocean  waves.  This  simplified 
hypothesis  ignores  all  the  complications  of  microseismic  generation,  prop¬ 
agation,  and  sensing.  However,  the  application  of  the  linear  assumption 
appeared  to  be  somewhat  successful  in  estimating  nearshore  waves  along 
the  Oregon  coast  as  reported  by  Creech  (1981)  and  Thompson,  Howell  and 
Smith  (1985). 

The  OSU  wave  meter  applies  a  linear,  u>~3  ,  transfer  function  to  the 
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seismometer  velocity  signal  and  displays  the  resulting  voltage  as  a  wave 
amplitude.  The  voltage  is  empirically  calibrated  with  a  square  root  law, 
arising  from  the  Longuet-Higgins’  generation  theory.  Ocean  wave  periods 
are  assumed  to  be  exactly  double  the  observed  periods  of  the  wave  meter 
output. 

In  the  present  work,  we  also  investigated  two  potential  improvements 
to  the  simple  assumptions  of  the  wave  meter.  First,  the  use  of  spectral 
analysis  of  the  seismometer  signal  to  improve  estimates  of  the  dominant 
period  of  ocean  waves  was  considered.  Second,  the  use  of  the  primary 
inicroseismic  spectrum  was  investigated  to  assist,  in  correctly  estimating 
the  dominant  period,  and  to  determine  if  estimates  of  nearshore  waves 
could  be  separated  from  offshore  waves.  Since  the  primary  interest  of  the 
Corps  is  for  sea  conditions  in  coastal  areas  rather  than  in  deep  water,  the 
later  effort  was  considered  important.  Nearshore  wave  characteristics  often 
vary  considerably  from  those  in  deep  water  because  of  the  complexities 
of  shallow- water  wave  propagation,  refraction,  diffraction,  reflection,  and 
shoaling.  The  local  bathymetry,  wind  and  currents,  as  well  as  the  frequency 
and  directional  spectrum  of  the  offshore  waves  affect  the  variations.  For 
this  reason,  a  direct  measurement  of  the  shallow-water  waves  is  usually 
required. 

The  original  design  of  the  wave  meter  relies  on  manual  strip  chart  anal¬ 
ysis  to  determine  wave  periods.  An  obvious  improvement  would  be  the 
digitizing  output  signals  from  the  wave  meter  and  the  analyzing  time-series 
periods  by  a  computer.  This  approach  has  been  implemented  on  an  exper¬ 
imental  basis  as  part  of  the  work  reporter!  here.  The  period  was  estimated 
using  an  average  zero- crossing  algorithm  to  approximate  the  manual  anal¬ 
ysis  method.  This  time-domain  analysis,  however,  can  only  he  expected  to 
yield  period  estimates  comparable  to  measurements  from  wave  spectra,  if 
the  spectral  shape  of  the  seismometer  signal  is  similar.  Longuet-Higgins" 
analysis  (1950)  was  based  on  monochromatic,  deterministic  waves  and  not 
a  random  sea  surface.  Hasselmanri  (19C3)  derived  several  cases  of  spectral 
shape  of  the  seismic  response,  which  were  not,  in  general,  the  same  as  the 
wave  spectrum.  Kadota  and  Labianca  (1981)  proposed  that,  in  deep  wa¬ 
ter.  the  spectral  shape  of  the  seismic  response  is  dependent  only  upon  the 
swell  portion  of  the  surface  wave  spectrum,  and  the  random  phase,  wind- 
sea  makes  no  contribution.  Their  formulation  predicts  a  spectral  shape 
of  the  seismic  reponse  which  is  different,  from  the  sea-surface  spectrum. 
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Darbyshiio  (1950)  and  Haubrich,  Munk  and  Snodgrass  (1963)  presented 
seismic  spectra  having  a  peak  at  the  frequency  which  is  twice  that  of  dom¬ 
inant  ocean  wave  frequency.  Especially,  Haubrich,  Munk  and  Snodgrass 
(1963)  pointed  out  the  correlation  between  the  dispersive  changes  of  the 
seismic  and  ocean  wave  peak  frequencies  with  time. 

Therefore,  it  is  reasonable  to  propose  that  a  seismic  wave  meter  system 
should  use  spectral  analysis  to  determine  period  estimates  rather  than  a 
zero-crossing  analysis.  In  the  present  work,  the  use  of  the  peak  of  the 
seismic  spectrum  is  examined  to  estimate  the  spectral  peak  of  ocean  weaves. 
Following  the  deterministic  analysis  of  Longuet-Higgins,  the  peak  of  the 
wave  spectrum  should  be  one  half  the  frequency  of  the  peak  of  the  seismic 
spectrum.  Under  the  more  general  assumptions  of  Hasselmann  (1963)  and 
Kadota  and  Labianca  (1981).  the  ratio  of  spectral  peaks  could  vary.  For 
a  swell-dominated  sea,  such  as  often  characterizes  the  high  sea  states  near 
the  Oregon  coast,  the  ratio  should  closely  approximate  two. 

The  "primary”  spectrum  of  microseisms  refers  to  small  amounts  of  en¬ 
ergy  which  appear  at  the  same  frequency  as  nearshore  ocean  waves.  This 
sort  of  longer  period  seismic  wave  is  also  well  known  from  experiments,  for 
example,  Pomeroy  ( 1959)  and  Oliver  (1962).  These  seismic  waves  appear  to 
have  the  same  period  as  ocean  waves.  The  spectral  results  of  microse'sms 
presented  by  Haubrich,  Munk  and  Snodgrass  (1963)  showed  that  the  micro¬ 
seisms  recorded  near  San  Diego,  California,  had  a  primary  frequency  peak, 
of  small  intensity,  but  clearly  related  to  ocean  swell.  Thus,  these  earlier 
studies  accordingly  suggested  a  theory  that  such  seismic  activities  must  be 
due  to  local  water  waves  impacting  a  ''oastline.  This  theory  was  extended 
to  an  extensive  mathematical  development  by  Hasselmann  (1963). 

The  existence  of  the  primary  spectrum  related  to  nearshore  waves  could 
provide  a  mechanism  to  solve  one  of  the  problems  of  the  existing  measure¬ 
ment  system  of  seismic  waves.  That  is  the  inability  to  separate  seismic 
responses  to  ocean  waves  of  interest  locally  from  those  to  other  sources, 
for  example,  ocean  waves  at  a  great  distance  from  the  seismometer.  A  sec¬ 
ondary  application  could  be  to  provide  assistance  in  estimating  the  period 
of  the  local  sea  state. 
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Part  2 

Experimental  Data  Collection 


2.1  Experimental  Plan 

The  experimental  plan  was  conceived  as  a  minimum  cost  effort  to  acquire 
data  sufficient  to  examine  the  proposed  improvements  to  the  seismic-wave 
measurement  system.  An  existing  OSU  seismic-wave  meter  system  at  the 
Chetco  River  Coast  Guard  Station  near  Brookings,  Oregon,  was  retrofitted 
with  a  digital  data  logger.  The  data  acquisition  period  was  planned  to 
coincide  with  wave  measurements  by  other  groups  in  order  to  provide  com¬ 
parison  data.  Figure  1  shows  the  locations  of  the  various  measurement 
sites. 


2.2  Seismic  Wave  Meter  Data  Acquisition 

The  seismic  wave  meter  is  composed  of  a  Teledyne- Geotech  Model  SL-210 
portable  seismometer  designed  for  geophysical  surveys.  As  used  with  the 
wave  meter,  it  is  adjusted  for  critical  damping  with  a  natural  period  of 
about  18  sec.  The  output  of  the  seismometer  is  connected  to  a  preamplifier 
and  filter  which  produces  the  wave  meter  output.  For  this  experiment  the 
output  of  the  seismometer  and  the  filter  were  simultaneously  recorded  by 
a  Sea  Data  Model  1250  digital  data  logger.  The  signals  were  tapped  from 
the  wave  meter  circuitry  which  remained  in  operation.  They  were  then 
amplified  and  scaled  to  within  the  range  of  the  data  logger  by  a  signal 
conditioner  composed  of  adjustable  instrumentation  amplifiers.  Figure  2  is 
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a  block  diagram  of  the  data  acquisition  system.  Data  from  both  channels 
were  sampled  at  2  Hz  for  17  min  every  4  hr. 

The  wave  meter  is  normally  empirically  calibrated  to  observed  wave  con¬ 
ditions  at  each  site.  The  calibration  for  the  wave  meter  voltage  recorded 
by  the  data  logger  was  determined  by  adjusting  the  wave  meter  zero  con¬ 
trol  over  the  operating  range  of  the  meter  indicator,  while  simultaneously 
measuring  the  voltage  input  to  the  data  logger  and  manually  recording 
the  wave-height  indication  on  the  meter.  Since  the  meter  scale  is  used  to 
provide  the  square  root  dependency  of  the  output,  the  resulting  curve  is  a 
bipolar  quadratic.  Figure  3  is  a  plot  of  a  curve  fit  through  the  measured 
values  and  is  the  calibration  used  for  the  wave-meter  data  analysis.  It  could 
not  be  established  that  the  wave  meter  site  at  Chetco  River  had  ever  been 
formally  calibrated  to  observed  wave  conditions.  No  recalibration  or  ad¬ 
justment  of  the  “K”  value  had  been  performed  in  the  recent  memory  of 
personnel  familiar  with  the  installation.  The  seismometer  and  electronics, 
however,  had  been  regularly  maintained  and  were  in  good  operating  condi¬ 
tion  both  during  and  after  the  experiment.  Because  of  the  question  about 
the  absolute  calibration,  interpretation  of  the  magnitudes  of  comparison 
data  should  be  done  in  a  relative  sense. 

No  attempt  was  made  to  provide  an  absolute  calibration  of  the  raw 
seismic  signal.  The  velocity  signal  was  amplified  by  an  arbitrary  gain  to 
provide  sufficient  sensitivity  to  microseisms.  The  amplified  signal  was  fil¬ 
tered  by  a  low  pass,  6  pole,  Butterworth  filter  with  a  corner  frequency  of  1 
Hz.  The  resulting  voltage  was  recorded  by  the  data  logger. 

2.3  Comparison  Data 

For  comparison,  wave  data  of  shallow  water  were  provided  from  an  acoustic 
current  meter  located  in  78  ft  of  water  just  offshore  from  the  seismometer 
site.  This  gage  was  part  of  a  current  meter  deployment  by  Dr.  Charles  K. 
Sollitt,  OSU,  under  the  sponsorship  of  the  Portland  District. 

Since  data  from  OSU  were  not  always  available  during  comparison  peri¬ 
ods,  supplementary  comparison  sites  were  also  used.  Near  Coquille  River, 
a  waverider  buoy  in  deep  water  and  an  Sxy  array  in  shallow  water  are  op¬ 
erated  for  the  Corps,  Field  Wave  Gaging  Program  by  the  Scripps  Institute 
of  Oceanography  (SIO).  Although  these  gages  are  located  some  distance 
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from  the  Chetco  seismometer,  they  are  exposed  to  similar  offshore  wave 
conditions  (Figure  1). 
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Part  3 

Data  Analysis  and 
Comparisons 


3.1  Evaluation  of  Estimates  of  Significant  Wave 
Height  by  Existing  Wave  Meter 

Figures  4  through  9  show  the  time-series  of  significant  wave  heights  from 
the  seismometer  and  comparison  gages.  The  significant  wave  heights  fur¬ 
nished  by  SIO  and  OSU  were  calculated  from  spectral  estimates  as  four 
times  the  square  root  of  zero  moment.  The  significant  wave  heights  for  the 
seismometer  were  calculated  by  averaging  the  highest  one-third  waves  in 
the  time  domain. 

In  Figure  4,  the  periods  from  Julian  day  (J.D.)  195  through  198  show 
comparable  estimates  from  both  the  seismometer  and  the  closely  located 
OSU  gage.  From  J.D.  198  through  214  (Figure  5),  however,  the  estimates 
begin  to  diverge  widely.  In  fact,  the  seismometer  data  show  better  agree¬ 
ment  with  the  deepwater  buoy  data  of  SIO  during  this  period  than  the 
closer,  but  shallower  data  of  OSU. 

Again,  during  two  moderate  storms  in  August  (Figure  6),  it  can  be  seen 
that  the  seismometer  estimates  follow  the  general  pattern  of  the  deepwa¬ 
ter  buoy.  Unfortunately,  Chetco  River  data  were  not  available  during  this 
time  period.  While  a  systematic  difference  of  the  estimates  would  not  be 
surprising  given  the  questions  about  calibration  of  the  seismometer  gage, 
the  seismometer  appears  to  overestimate  the  wave  heights  during  the  first 
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storm  and  underestimate  during  the  smaller  second  storm.  Without  com¬ 
parison  data,  it  cannot  be  ruled  out  that  the  area  near  the  seismometer 
responded  differently  to  the  two  storms. 

Figures  7  and  8  show  results  from  two  winter  storms.  Again  Chetco 
River  data  were  not  available,  but  good  agreement  can  be  seen  in  the  shape 
of  the  variation  of  significant  wave  heights  during  the  course  of  the  storms, 
except  for  a  systematic  difference  in  the  amplitude. 


3.2  Evaluation  of  Period  Estimates  by  Ex¬ 
isting  Wave  Meter 

The  results  for  time-series  of  wave  periods  are  presented  in  Figures  10 
through  15.  Wave  periods  for  OSU  and  SIO  estimates  represent  the  period 
at  which  a  maximum  density  occurs  in  the  wave  spectrum.  For  the  wave- 
meter  estimates,  the  average  zero-crossing  period  of  the  highest  one-third 
waves  was  calculated  to  be  consistent  with  previous  manual  methods  used. 
The  seismometer  data  are  plotted  without  applying  the  constant  multiplier 
of  two.  The  wave  periods  obtained  from  OSU  and  SIO  may  be  used  to 
estimate  the  significant  wave  periods  by  using  one  of  the  experimental  for¬ 
mulae  (see  Goda  1974).  However,  this  would  result  in  only  slight  changes 
in  their  values. 

These  plots  show  comparisons  of  the  period  estimates  from  the  same 
data  sets  as  the  significant  wave-height  comparisons.  Based  on  traditional 
wave  meter  analysis  procedures,  the  wave  meter  period  should  be  one  half 
the  wave  period. 

The  periods  of  typical  low  summer  waves  such  as  those  in  Figures  10 
and  11  show  some  events  where  the  wave  meter  estimates  are  reasonably 
close  to  the  periods  measured  by  the  OSU  gage.  There  are,  however,  long 
sequences  where  the  seismometer  wave  meter  appears  to  be  responding  to  a 
period  quite  different  from  those  of  the  OSU  gage.  Unlike  the  wave-height 
comparisons,  the  period  estimates  do  not  show  any  strong  correlation  with 
the  SIO  deepwater  buoy.  During  the  storm  events  shown  in  Figures  12 
through  14,  period  estimates  come  much  closer  to  one  half  of  the  observed 
wave  periods.  Although  the  apparent  correlation  is  much  stronger,  there 
are  still  several  segments  when  the  estimates  diverge  sharply.  Again,  since 


10 


shallow-water  data  near  the  seismometer  were  not  available,  local  effects 
could  not  be  evaluated. 


3.3  Spectral  Analysis  of  Raw  Seismometer 
Data 

Figures  16  through  20  show  smoothed  spectra  of  the  raw  seismometer  ve¬ 
locity  signal.  The  three-dimensional  plots  are  composed  of  a  sequence  of 
autospectra,  one  for  each  17-min  data  record.  The  raw  seismic  data  include 
the  same  measurement  periods  as  the  data  from  the  wave  meter  presented 
in  the  last  section. 

The  plots  from  the  July  1985  measurements  show  dominant  peaks  at  the 
secondary  frequencies.  The  shape  of  the  spectra,  however,  is  quite  broad 
and  contains  multiple  peaks.  Because  of  the  difficulty  in  applying  peak 
detection  algorithms  to  spectra  with  such  a  high  variation,  peak  frequency 
estimates  were  computed  at  the  centroid  of  the  spectra,  or  the  period  of 
the  first  moment.1 

It  is  clear  from  the  plots  that  an  energetic  secondary  spectrum  exists 
for  every  wave  record,  both  for  the  small  wave  heights  during  the  summer 
conditions,  and  for  the  storm  events  in  August  1985  and  March  and  April 
1986  (Figures  18  through  20).  A  microcomputer  analysis  program  run¬ 
ning  in  real-time  should  be  able  to  robustly  compute  spectra  from  a  raw 
seismometer  signal  and  compute  period  estimates  from  the  centroid  of  the 
spectrum. 

Figures  21  and  22  show  the  results  of  period  calculations  from  the  spec¬ 
tra  together  with  SIO  estimates  for  the  August  storm.  The  general  evo¬ 
lution  of  the  seismic  period  estimates  follows  closely  the  measured  data, 
especially  the  deepwater  buoy.  The  variations  in  period  show  the  expected 
dispersive  effects  of  storm  wave  arrival  which  has  been  noted  in  previous 
seismic  investigations  (Haubrich,  Munk  and  Snodgrass  (1963)). 

Figures  23  and  24  show  a  scatter  diagram  of  periods  estimated  from 


7  IfS(f)df 
1  fS(f)df  ' 

where  f  =  peak  frequency  estimate,  /  =  frequency,  and  S(f )  =  power  spectral  density. 
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the  seismic  spectra  versus  the  comparison  periods  for  both  the  deepwater 
data  and  shallow-water  data  from  SIO.  The  plot  includes  all  data  points 
from  both  the  winter  and  summer  cases.  The  grouping  of  the  points  is 
quite  good,  especially  with  the  deepwater  data.  The  slope  of  the  linear 
regression  line  through  the  points  is  0.47.  The  scatter  of  points  with  the 
shallow  water  comparison  data  is  not  as  tight,  and  the  slope  of  the  linear 
regression  line  is  0.42.  Again,  the  distance  of  this  shallow-water  data  from 
the  seismometer  site  could  account  for  these  differences. 

To  examine  further  the  implications  of  using  a  simple  double-frequency 
assumption  to  estimate  wave  period  from  the  seismometer  spectrum,  in 
Figures  25  and  26  we  plot  a  histogram  of  the  ratio  of  comparison  data  to 
the  seismometer  spectral  peak  period.  For  the  deepwater  data,  the  mean 
value  is  2.16  with  a  standard  deviation  of  0.27.  When  the  ratio  is  computed 
with  the  shallow-water  data,  the  mean  is  2.37,  and  the  standard  deviation 
is  0.3. 

3.4  Investigation  of  Primary  Frequency  Com¬ 
ponents  in  Seismic  Spectra 

Close  inspection  of  the  microseismic  spectral  plots  for  the  summer  data 
(Figures  16  through  18)  show  clear  evidence  of  small  amounts  of  energy  at 
frequencies  lower  than  the  secondary  frequency  peaks.  The  energy  appears 
around  0.15  Hz,  or  close  to  the  6-  to  7-sec  periods  reported  from  the  SIO 
deepwater  buoy  as  shown  previously  in  Figures  10  and  11.  Since  this  fre¬ 
quency  is  also  approximately  one  half  the  peak  periods  estimated  from  the 
secondary  frequencies,  it  may  be  concluded  that  these  peaks  are  primary 
frequency  energy.  It  is  also  interesting  to  observe  that  the  relative  amount 
of  energy  in  the  primary  peaks  does  not  always  correlate  with  the  energy 
in  the  secondary  peaks,  supporting  the  hypothesis  of  previous  investigators 
that  different  physical  generation  processes  are  involved  in  the  primary 
and  secondary  energy.  The  fact  that  primary  spectral  energy  appears  to  be 
present  much  of  the  time,  even  in  relatively  mild,  summer  sea  conditions  is 
also  an  important  observation. 

Figures  27  and  28  are  expanded  portions  of  the  spectra  from  part  of 
the  July  data.  At  this  scale,  energy  in  the  period  range  around  0.06  Hz, 
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which  appeared  flat  in  the  previous  plots,  is  clearly  visible.  Returning  to 
the  comparison  between  seismometer  period  estimates  and  the  OSU  gage 
near  Chetco  River  shown  in  Figures  10  and  11,  there  was  a  time  starting 
about  July  16  where  the  seas  measured  by  the  pressure  gage  appeared  to 
come  under  the  influence  of  a  swell  of  about  17  sec.  During  this  period 
the  estimates  from  the  wave  meter  became  somewhat  confused.  It  does 
appear  in  comparing  Figures  27  and  10  that  observable  energy  at  the  16- 
to  17-sec  period  in  the  seismic  spectrum  correlates  with  the  times  that  the 
pressure  gage  was  showing  similar  periods.  Also  during  this  period,  the 
SIO  deepwater  buoy  continued  to  show  the  periods  of  the  wind-sea. 

Figure  29  is  a  single  spectrum  of  the  seismometer  signal  from  July  19, 
when  the  swell  part  of  the  primary  spectrum  was  particularly  energetic. 
This  plot  shows  some  of  the  opportunities  as  well  as  the  potential  problems 
in  considering  the  primary  spectrum.  Since  the  sea  state  was  composed 
of  both  a  swell  component  at  17  sec  and  a  wind-sea  at  about  7  sec,  these 
components  can  be  identified  as  having  energy  in  both  the  primary  and 
secondary  frequencies.  Unfortunately,  the  secondary  swell  energy,  theoret¬ 
ically  falling  at  0.12  Hz,  comes  very  close  to  the  primary  wind-sea  at  0.14 
Hz.  This  could  present  problems  in  separating  the  contributions  without 
prior  knowledge  of  the  spectrum. 

Turning  to  the  case  of  winter  data,  the  spectral  densities  of  the  sec¬ 
ondary  frequencies  are  more  than  an  order  of  magnitude  larger  than  the 
the  summer  spectra.  In  most  cases  any  primary  energy  is  completely  ob¬ 
scured  at  the  plotting  scale  of  Figures  18  through  20.  For  example,  Figure 
30  shows  a  single  spectrum  plot  from  March  11  in  a  form  similar  to  Figure 
29.  No  primary  energy  is  visible.  Looking  at  the  August  storm,  Figure  31 
shows  an  expanded  scale  of  the  spectrum  from  August  2.  Again  very  little 
primary  energy  is  visible.  Eight  hours  later,  however,  Figure  32  shows  a 
peak  has  developed  around  0.07  Hz.  This  corresponds  with  the  develop¬ 
ment  of  the  storm  and  the  gradual  shift  towards  a  swell  period  as  measured 
by  the  SIO  buoy  as  shown  in  Figures  6  and  12.  The  same  development  can 
be  seen  for  the  March  storm  in  Figures  33  and  34.  The  shape  of  the  primary 
spectrum  in  Figure  34  looks  remarkably  similar  to  a  wave  spectrum. 
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3.5  Wave-Height  Estimates  from  Microseis- 
mic  Spectra 

Spectral  analysis  of  the  raw  seismic  signal  would  permit  estimation  of  sig¬ 
nificant  wave  heights  using  the  same  assumptions  and  the  equivalent  signal 
processing  as  employed  by  the  linearizer  circuit  of  the  wave  meter.  A  digital 
low  pass  filter  with  the  same  response  characteristics  could  be  applied  and 
the  results  either  analyzed  in  the  time  domain  as  in  the  evaluation  section 
of  this  report,  or  equivalent  estimates  could  be  computed  based  on  the  zero 
moment  of  the  filtered  microseismic  spectra.  In  any  case,  the  results  would 
be  the  same  as  previously  presented. 

Potential  improvements  to  wave-height  estimates  could  be  pursued  by 
attempting  to  evaluate  the  energy  of  the  primary  and  secondary  frequen¬ 
cies.  It  is  hoped  that  the  primary  frequencies  may  assist  in  separating  the 
effects  of  nearshore  and  distant  seismic  activity.  Some  encouragement  was 
found  for  the  case  of  periods  above.  Unfortunately,  Loo  little  nearshore 
comparison  data  were  available  to  evaluate  use  of  the  primary  frequency 
energies.  Simultaneous  offshore,  nearshore,  and  seismometer  data  for  a  long 
measurement  period  would  be  required.  Given  that  adequate  correlation 
between  primary  peak  energy  and  local  significant  wave  height  could  be 
established,  applying  this  in  a  wave  measurement  system  would  present 
additional  difficulties.  Because  of  the  very  small  energies  of  the  secondary 
peaks,  and  the  possible  occurrences  of  overlap  as  noted  above,  algorithms 
for  separation  of  the  spectra  would  be  more  difficult  than  the  relatively 
straightforward  computations  used  for  the  period  estimates.  The  acqui¬ 
sition  of  the  needed  data  and  the  development  of  such  algorithms  could 
present  a  useful  path  for  future  work. 

For  a  simple  analysis  procedure  based  on  the  microseismic  spectrum, 
the  following  was  considered.  Remembering  that  the  primary  spectrum 
contributes  very  little  energy  to  the  total,  some  simple  comparisons  were 
computed.  The  total  energy  of  the  raw  seismic  spectrum  was  computed  as 
the  square  root  of  the  zero  moment.  Figures  35  and  36  show  the  square  root 
of  the  total  energy  plotted  against  observed  wave  height  for  both  deepwater 
and  shallow-water  data  from  SIO.  For  the  deep-water  scatter  diagram,  the 
points  fall  reasonably  close  to  a  straight  line,  except  for  a  number  of  records 
when  the  wave  height  was  grossly  under  predicted.  Because  of  the  limited 
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amount  of  comparison  data  close  to  the  seismometer,  it  is  difficult  to  make 
conclusions  based  on  these  results,  but  the  agreement  does  not  seem  to  be 
significantly  worse  than  that  obtained  from  the  present  wave  meter  method. 
Of  course,  use  of  either  of  these  methods  would  still  rely  on  the  site-specific 
calibration  of  the  curve  as  is  presently  the  case  for  the  wave  meter. 
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Part  4 

Conclusions  and 
Recommendations 

4.1  Wave  Meter  Accuracy 

The  results  of  comparing  significant  wave  height  estimates  from  the  existing 
wave  meter  system  for  summer  wave  conditions  were  not  as  favorable  as 
the  results  previously  reported  for  winter  storms  only.  Results  from  the 
small  amount  of  winter  data  compared  here  show  greater  correlation  than 
from  the  summer  data,  but  large  errors  existed  in  actual  estimates.  Better 
site-specific  calibration  of  the  seismometer  at  Chetco  River  could  reduce 
these  errors.  As  in  the  previous  report,  period  estimates  frequently  had 
large  errors. 

A  limited  amount  of  comparison  data  was  available  for  local  wave  con¬ 
ditions  near  Chetco  River.  No  local  data  were  available  for  the  winter 
comparison  periods.  Based  on  this  limited  set  of  summer  data,  the  wave 
meter  estimates  appear  to  be  more  strongly  influenced  by  distant  deepwater 
conditions  than  local  conditions.  For  cases  of  large  swell-dominated  storms, 
local  conditions  may  closely  parallel  distant  conditions. 

The  seismic  system  operated  reliably  during  the  entire  data  acquisition 
period  and  provided  continuous  data  except  during  periods  of  geoseismic 
activity.  Periods  of  seismic  activity  which  prevent  wave  estimates  can  last 
from  hours  to  days. 
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4.2  Improvements  to  the  Seismic  Wave  Mea¬ 
surement  System 

Directly  digitizing  and  processing  the  raw  seismic  signal  are  feasible.  Spec¬ 
tral  analysis  of  the  seismic  data  improves  estimates  of  wave  period.  Initial 
approaches  should  use  the  energetic  secondary  frequencies  to  estimate  the 
peak  period  and  then  use  a  constant  multiplier  of  two  to  estimate  the  ocean 
wave  period.  Additional  comparison  data  and  investigation  may  allow  the 
use  of  the  primary  frequency  to  assist  in  separating  the  periods  at  the  local 
site  from  more  distant  conditions.  Estimates  of  significant  wave  height  can 
be  made  using  a  digital  system.  These  estimates  should  be  at  least  as  good 
as  the  existing  wave  meter. 

It  is  well  within  the  state  of  the  art  to  design  and  construct  a  ruggedized 
microcomputer  system  capable  of  digitizing  and  processing  the  seismometer 
signal  in  real-time.  Estimates  of  wave  height  and  period,  along  with  qual¬ 
ity  control  checks,  could  be  produced  using  analysis  methods  similar  to 
those  reported  here.  Wave  estimates  could  be  stored  internally  and  trans¬ 
mitted  to  a  central  site  via  telephone,  radio,  or  satellite  communications 
technology.  The  microcomputer  system  could  be  constructed  to  withstand 
a  severe  marine  environment.  The  seismometer  itself,  however,  requires  a 
temperature-controlled  environment.  The  microcomputer  system  could  be 
designed  such  that  its  power  consumption  would  be  low  enough  to  allow 
solar-powered  operation.  The  power  required  to  temperature  control  the 
seismometer  would  be  the  limiting  factor  on  total  system  power  consump¬ 
tion.  In  most  cases,  the  seismometer  heating  requirements  would  require 
connection  to  an  electric  utility  or  some  other  source  of  accurately  con¬ 
trolled  heat.  The  development  of  real-time,  robust  versions  of  the  analysis 
programs  would  represent  the  major  development  cost  of  the  system. 
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PACIFIC  OCEAN 


Figure  1.  Measurement  site  locations 


Figure  2.  Seismometer  data  acquisition 
block  diagram 
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Figure  3.  Calibration  curve  for  seismic  data  from 
data  logger  voltage  to  wave  height 
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Figure  4.  Time  series  of  significant  wave  heights 
from  11-21  July  1985 
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Figure  5.  Time  series  of  significant  wave  heights 
from  23  July  to  I  August  1985 
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Figure  6.  Time  series  of  significant  wave  heights 
from  2-12  August  1985 
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Figure  10.  Time  series  of  significant  wave  periods 
from  11-21  July  1985 


Figure  11.  Time  series  of  significant  wave  periods 
from  22  July  to  1  August  1985 
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Figure  13.  Time  series  of  significant  wave  periods 
from  9-10  January  1986 
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Figure  14.  Time  series  of  significant  wave  periods 
from  10-14  March  1986 
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Figure  20.  Frequency  spectra  of  microseisms 
from  11-14  April  1986 
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Figure  22.  Dispersive  change  of  peak  frequency  in  seismic  and 
ocean  wave  spectra  at  water  depth  of  52  ft 
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Figure  25.  Microseisms  frequency  versus  ocean  wave  frequency.  Ratio  of  peak  spectrum 
frequencies  at  water  depth  of  230  ft  from  July  and  August  1985  and  March  and  April  1986 
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Figure  26.  Microseisms  frequency  versus  ocean  wave  frequency.  Ratio  of  peak  spectrum 
frequencies  at  water  depth  of  230  ft  from  July  and  August  1985  and  March  and  April  1986 
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Figure  29.  Frequency  spectrum  of  microseisms 
from  19  July  1985  at  0630 
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Figure  30.  Frequency  spectrum  of  microseisms 
from  11  March  1986  at  1007 
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Figure  31.  Frequency  spectrum  of  microseisms 
from  2  August  1985  at  1430 
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Figure  32.  Frequency  spectrum  of  microseisms 
from  2  August  1985  at  2230 


0  08 


o©s/(2#*s;i°a)  'A^suep  jenoeds 


Frequency,  1/sec 

Figure  33.  Frequency  spectrum  of  microseisms 
from  10  March  1986  at  1807 
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Figure  34.  Frequency  spectrum  of  microseisms 
from  10  March  1986  at  2207 
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Figure  35.  Microseisms  versus  ocean  wave  heights  at  water  depth  of  230  ft 
from  July  and  August  1985  and  March  and  April  1986 
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Figure  36.  Microseisms  versus  ocean  wave  heights  at  water  depth  of  52  ft 
from  July  and  August  1985  and  March  and  April  1986 


Appendix  A 

CORRELATION  BETWEEN 
WAVE  HEIGHTS  AND 
ESTIMATES  BY 
SEISMOMETER 

Figure  A-l  shows  the  relationship  between  wave  height  estimates  from  the 
seismometer  and  significant  wave  heights  at  Coquille  River  (offshore  deep 
water)  from  SIO.  The  relationship  between  wave-height  estimates  from  the 
seismometer  and  significant  wave  heights  at  Chetco  River  from  OSU  is  not 
shown  here  because  the  correlation  is  very  poor  as  can  be  seen  in  figures  4 
and  5. 


